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Abstract 

We use a fractional transformation to connect the traveling wave solutions of the nonlinear 
Schrodinger equation (NLSE), phase-locked with a source, to the elliptic functions satisfying, /" ± 
a/± A/ 3 = 0. The solutions are necessarily of the rational form, containing both trigonometric and 
hyperbolic types as special cases. Bright, and dark solitons, as also singular solitons, are obtained 
in suitable range of parameter values. 
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Much attention has been paid to the study of the externally driven NLSE, after the 
seminal work of Kaup and Newell jlj. This equation features prominently in the problem 
of optical pulse propagation in asymmetric, twin-core optical fibers @, H HI ; currently an 
area of active research. Of the several applications of an externally driven NLSE, perhaps 
the most important ones are to long Josephson junctions charge densitv waves p, and 
the plasmas driven by rf fields Q. The phenomenon of autoresonance H, indicating a 
continuous phase locking between the solutions of NLSE and the driving field has been 
found to be a key characteristic of this system. In the presence of damping, this dynamical 
system exhibits rich structure including bifurcation. This is evident from analyses, around 
a constant background, as well as numerical investigations 0,0, U|. Although the NLSE 
is a well-studied integrable system no exact solutions have so far been found for the 



NLSE with a source, to the best of the authors' knowledge. All the above inferences have 
been drawn through perturbations around solitons and numerical techniques. 

In this Letter, we map exactly, the traveling wave solutions of the NLSE phase-locked with 
a source, to the elliptic functions, through a fractional transformation (FT). It was found that 
the solutions are necessarily of the rational type, with both the numerator and denominator 
containing terms quadratic in elliptic functions, in addition to having constant terms. It is 
well-known that the solitary wave solutions of the NLSE ["Lsl . 14l | are cnoidal waves, which 



contain the localized soliton solutions in the limit, when the modulus parameter equals one 
[Tif . Hence, the solutions found here, for the NLSE with a source, are nonperturbative in 
nature. We find both bright and dark solitons as also singular ones. Solitons and solitary 
pulses show distinct behavior. In the case, when the source and the solutions are not phase 
matched, perturbation around these solutions may provide a better starting point. 

For nonlinear equations, a number of transformations are well-known in the literature, 
which map the solutions of a given equation to the other 0, 17]. The familiar example 



being the Miura transformation |18j |. which maps the solutions of the modified KdV to those 
of the KdV equation. To find static and propagating solutions, appropriate transformations 
have also been cleverly employed, to connect the nonlinear equations to the ones satisfied 
by the elliptic functions: /" ± af ± A/ 3 = 0. Here and henceforth, prime denotes derivative 
with respect to the argument of the function. Solitons and solitary wave solutions of KdV, 
NLSE, and sine-Gordon etc., can be easily obtained in terms of the elliptic functions in this 
manner. 

We consider the NLSE with an external traveling wave source: 

i^ + ^+9\4\'4 + rt = kjM-»* , (1) 

where g, p, and k are real and £ = a(x — vt). The traveling wave solution, phase-locked with 
the external source is taken to be ip(x, t) = e 1 '*^ - "*] p(£) . Separating the real and imaginary 
parts of Eq. (1), and integrating the imaginary part one gets, 

where c is the integration constant. In order that, the external phase is independent of i/j, 
we put c = to obtain, 

a 2 p" + gp z - ep - k = 0, (3) 

2 

where e = u — \ — p. For illustrating our procedure, we assume that e > 0, and g > 0, 
the other cases can be studied analogously. After a suitable scaling of the field variable, the 
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above equation can be cast into the convenient form, 



p" + P 3 - P - k = 0. (4) 
We find, after a straightforward but lengthy algebra that, the following FT, 

for AS - A / 0, maps the solutions of Eq. (4) to the elliptic functions /(£), provided 
5 = 2. /(0 satisfies /" = / - / 3 , with a conserved quantity £ = f 2 /2 + (1/4)/ 4 - / 2 /2. 
The above equation has four bounded periodic solutions which are Jacobi elliptic functions: 
cn(£, m), sd(£,m), dn(£,m), and nd(£,m), where m is the modulus parameter. Apart from 
finite energy oscillatory solutions, the above has localized soliton solutions. It should be 
noted that for the attractive case (g < 0), the bounded solutions are sn(£, m), and cd(£,m). 
For AD = B, one only finds a constant amplitude solution. 

Before elaborating on specific number of interesting features emerging from the 

above mapping is worth mentioning. First of all, the nontrivial solutions are necessarily of 
the rational type, i.e., D ^ 0, in the presence of the source. Secondly, E = 0, and E ^ 
cases, show characteristically different behavior. For example, for Eq ^ case, one can have 
solutions with A = 0, and B ^ 0, of the type 

P(0 = 1 + (1/8 W <6> 

However, A ^ 0, and B = is not allowed. For E = case, one can have for B = 0, and 
A 7^ 0, a singular solution of the following form, 

P(0 = ^ (7) 

As we will see explicitly below, nonsingular solutions are also possible. Like the former case, 
here B^0, and A = is not allowed in the presence of the source. 

For explicitness, we consider the unsealed equation containing all the parameters and 
illustrate below various type of solutions, taking / = cn(£, m). Other cases can be similarly 
worked out. The consistency conditions are given by, 

Ae - 2a 2 (AD - B)(l - m) + gA 3 - k =0, (8) 
2eAD + eB + 6a 2 (AD - B)D(l - m) - 

Aa 2 (AD — B)(2m — 1) + 3gA 2 B — 3kD = 0, (9) 

AeD 2 + 2eBD + Aa 2 (AD - B)D(2m - 1) + 

6a 2 (AD - B)m + 3gAB 2 - 3kD 2 = 0, (10) 

eBD 2 — 2a 2 (AD — B)Dm + gB 3 — kD 3 = 0. (11) 

The above equations clearly indicate that the solutions, for m = 1, m = and other 
values of m, have distinct properties. For example, when m — 1, A is obtained as the 
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solution of the cubic equation [Eq. (8)], containing the source strength k. Similarly, for 
m — 0, either B, or D appears as the solution of Eq. (11). As noted earlier, when D — 0, 
B also equals zero, indicating only a constant solution (A ^ 0). One needs to be careful in 
choosing the real solutions of the above cubic equations, for suitable range of the parameter 
values. Although a wide class of solutions are allowed; for brevity, we only outline a few of 
the interesting solutions and their properties in a restricted range of equation parameters. 

We start with the localized solitons, by taking m = 1 (cn(£, 1) = sech(£)), with the 
parameter values A = 1, B = 5, and D = 1 — 5 we obtain 

l + feech 2 (Q 

P{0 = 1 + Fsech^)' (12) 

here the amplitude, width, and velocity are related as, 5 = —(Q + \JQ 2 — 4PR)/2Q, with 
Q = e - 2o? - 3k, P = 2e - 3k, R = -k - 2a 2 , and T = 1 + 5. Here the width a is the only 
independent parameter. As the above form of the solution indicates, both nonsingular and 
singular solitons are possible solutions depending on the values of e, and the source strength 
k. We have checked that periodic solutions are also allowed. Below we give a few more 
solutions, illustrating the differences between the localized soliton and the periodic solitary 
wave solutions. 

Case(I): Trigonometric solution.- In the limit m — 0, unlike the unperturbed NLSE, in 
this case, one finds rational solutions of the trigonometric type. Apart from the general 
solutions, interestingly, for these, one can obtain special cases, where A = and B ^ 
is allowed. However, the vice-versa is forbidden. Following is an example of this type of 
nonsingular periodic solutions, for the repulsive case: 

2k cos 2 (Q 

p(0 = ^i-fcos^r (13) 

Here, e has to be negative since, a 2 = — e/4; and is given by e = (—21gk 2 /2) 1 ^. This 
periodic solution is found to be stable, as evidenced from the numerical simulations seen in 
Fig.l. 

Case(II): Hyperbolic solution.- Unlike the above periodic case, here one finds that, the 
solutions with B = 0, and A / are allowed, the vice-versa not being true. Hence, 
these localized solitons behave differently from the soliton trains. There are both singular 
and nonsingular solutions. For B = 0, and m — 1, we found that a 2 = e/4, and e = 
( — 27(?fc 2 /2) 1 / 3 . This yields, the singular hyperbolic solution, 

3k 1 

The singularity here may correspond to extreme increase of the field amplitude due to self- 
focussing, as is known for the other nonlinear systems. We have also found nonsingular 
solutions of the above type. 

Case(III): Pure cnoidal solutions.- In the general case, for < m < 1, from Eq. (11), 
one obtains the width a, in terms of the other parameters. The other three equations then 
yield A, B, and D; in the process one encounters a cubic equation whose real roots need 
only to be considered. This puts restriction on the solution parameters. Interestingly, the 
constraints imply an m independent condition on the parameters in the form 

5eAD 2 + eAD 3 + eBD + 5eBD 2 + 3gA 2 BD 
+3gADB 2 - QkD 2 - QkD 3 + 3gA 3 D 2 + 3gB 3 = . 
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FIG. 1: Plot depicting the evolution of the trigonometric solution, for various times. 



Below we list a few special cases. For A — 0, D — 1, and m = 5/8; we found that for 
a 2 = (2/7)e, and e = 7( — (?A: 2 /18) 1 ' /3 ; this corresponds to attractive case. Explicitly, the 
solution is given by, 

14fc cn 2 (£,m) 
p(0 = ^l + cn^m)- (15) 

For A = 0, and m = 1/2; it is found that a 2 = and e = (—27gk 2 )^, for which 

2V3/c cn 2 (^,m) 

" K) = ( — >i + >»({,».)■ (18) 

Since the localized solitons are usually robust, we have performed numerical simulations 
to check the stability of the solutions pertaining to Case(I), i.e., the trigonometric solution. 
It is worth pointing out that the numerical techniques based on the fast Fourier transform 
(FFT) are expensive as they require the FFT of the external source. Hence, we have used 



the Crank-Nicholson finite difference method [19] to solve the NLSE with a source, which 
is quite handy, and unconditionally stable. The initial conditions chosen from the exact 
solution are knitted on a lattice with a grid size dx = 0.005, and dt = 5.0 x 10~ 6 . The 
evolution of this solution, as depicted in Fig.l indicate that it is a stable solution. 

In conclusion, we have used a fractional transformation to connect the solutions of the 
phase-locked NLSE with the elliptic functions, in an exact manner. The solutions are nec- 
essarily of the rational type that contain solitons, solitary waves, as also singular ones. Our 
procedure is applicable, both for the attractive and repulsive cases. Because of their exact 
nature, these will provide a better starting point for the treatment of general externally 
driven NLSE. Considering the utility of this equation in fiber optics and other branches of 
physics, these solutions may find practical applications. 

We acknowledge profitable discussions with Dr. E. Alam regarding the algorithm that 
has been implemented here, and Prof. A. Khare for many useful discussions. 
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